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Abstract 

We investigate the geometric properties of the equilibrium manifold of a thermodynamic system 
determined by the van der Waals equations of state. We use the formalism of geometrothermody- 
namics to obtain results that are invariant under Legendre transformations, i. e., independent of 
the choice of thermodynamic potential. It is shown that the equilibrium manifold is curved with 
curvature singularities situated at those points where first order phase transitions occur. Moreover, 
the geodesic equations in the equilibrium manifold are investigated numerically to illustrate the 
equivalence between geodesic incompleteness and curvature singularities as a criterion to define 
phase transitions in an invariant manner. 
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I. INTRODUCTION 



It is well known that in nature there exist four different interactions and all of them 
can be interpreted in terms of concepts of differential geometry. Indeed, Einstein proposed 



the astonishing principle "field strengt 
gravitational field (see, for instance, |]J 



;h = curvature" to understand the physics of the 
1 

2|). In this case, curvature means the curvature of 
spacetime, a Riemannian manifold with a metric compatible with a connection that is unique 
as a consequence of the assumption that the torsion tensor vanishes. The second element of 
general relativity are Einstein's field equations R^ u — \g^ u R = 8nGT^ v that established for 
the first time the amazing principle "geometry = energy" . 

The conceptual fundamentals of this principle were very controversial when first for- 
mulated; however, experimental evidence has shown its correctness and even modern gen- 
eralizations of Einstein's theory follow the same principle. On the other hand, since the 
field strength can be considered as a measure of the gravitational interaction, we conclude 
that the entire idea of general relativity can be summarized in the principle "interaction = 
curvature" . 

The discovery by Yang and Mills [3| that the theory of electromagnetism can be inter- 
preted as an Abelian gauge theory and that non-Abelian generalizations can be constructed 
in a similar way, represents a major achievement. Today, it is known that electromagnetism 
can be described geometrically in terms of the elements of a principal fiber bundle. In this 
case, the base manifold is the Minkowski spacetime, the standard fiber is the gauge group 
U (1), which represents the internal symmetry of electromagnetism, and the connection across 
the fibers is a local cross-section that takes values in the algebra of U(l). 

This result opened the possibility of fixing the background metric, in accordance with the 
desired properties of the base manifold, and selecting different connections as local cross- 
sections of the principal fiber bundle. In particular, it was shown that connections with 
values in the Lie algebra of the gauge groups SU(2) and SU(3) can be used to represent 
the internal symmetries of the weak and strong interactions, respectively. This interesting 
geometrical approach constituted the basis for the construction of modern gauge theories 
that are used as the starting point to formulate the physics of the electromagnetic, weak, 
and strong interactions. It follows that the principle "curvature = interaction" holds for 
all known forces of nature, and is a fundamental element for the construction of the corre- 
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sponding gauge theories. 

Consider now the case of a thermodynamic system. In very broad terms, one can say that 
in a thermodynamic system the constituents are subject to the action of all known forces. 
However, it is not necessary to consider all the physical details of each interaction. In fact, 
due to the large number of particles involved in the system, only a statistical approach is 
possible from which average values for the physical quantities of interest, like energy, entropy, 
etc., are derived. The question arises whether it is also possible to find a geometric construc- 
tion for which the principle "curvature = thermodynamic interaction" holds. We will see 
in the present work that the formalism of geometrothermodynamics (GTD) 4] satisfies this 
requirement. First, we must mention that our interpretation of thermodynamic interaction 
is based upon the statistical approach to thermodynamics in which all the properties of 
the system can be derived from the explicit form of the corresponding Hamiltonian |5]. As 
usual, the interaction between the particles of the system is described by the potential part 
of the Hamiltonian. Consequently, if the potential vanishes, we say that the system does 
not possess thermodynamic interaction. 

In this work, we present the formalism of GTD using Riemannian contact geometry for 
the definition of the thermodynamical phase manifold and the equilibrium manifold. It is 
explained why it is necessary to introduce the auxiliary phase manifold in order to take into 
account the Legendre transformations in the appropriate manner. We apply our formalism 
to study the geometric properties of the equilibrium manifold of the van der Waals gas, and 
find that in general it is curved for arbitrary values of the parameters, except in the limiting 
case of an ideal gas for which the manifold turns out to be flat. Moreover, it is shown that 
the curvature singularities are located at those points where first order phase transitions 
occur. In addition, we explore the behavior of the geodesies, which represent quasi-static 
processes connecting different equilibrium states, and show that geodesic incompleteness can 
also be used as a criterion to detect points of phase transitions. 

This paper is organized as follows. In Section [III we introduce the main concepts of 
Riemannian contact geometry that are necessary to define the phase manifold and its equi- 
librium submanifolds. A particular metric is given that is invariant under total Legendre 
transformations, and is used to describe systems in which first order phase transitions can 
take place. Section DTD is dedicated to the investigation of the geometric properties of the 
equilibrium manifold for the van der Waals system, its curvature and singularities. Section 
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HVl includes an analysis of the geodesic equations that are integrated numerically to find the 
points where geodesic incompleteness appears. Finally, Section [V] is devoted to discussions 
of our results. Throughout this paper we use units in which k B = 1. 

II. BASIC ASPECTS OF GEOMETROTHERMODYNAMICS 

From a geometric point of view, the most important object of a given thermodynamic 
system is its thermodynamic space £, each point of which represents an equilibrium state. To 
be more specific it is convenient to introduce coordinates, say E a , on £ , where a — 1, 2, • • • , n. 
The integer n represents the number of thermodynamic degrees of freedom of the system. In 
general, it is possible to choose any set of n coordinates for £, but for the sake of convenience, 
the set {E a } of extensive thermodynamic variables is usually taken as the coordinates of the 
equilibrium space. Furthermore, a given thermodynamic system is uniquely characterized 
by the fundamental equation $ = Q(E a ), where $ is the thermodynamic potential that is 
used to describe the system Q]. The importance of the fundamental equation lies in the 
fact that from it one can derive all the equations of state that characterize the system in a 
unique manner. It follows that the fundamental equation can distinguish between different 
systems. 

Of course, for a fundamental equation to be physically relevant it must be in accordance 
with the laws of classical thermodynamics. For instance, the functional dependence of $ 
must be such that the second law is satisfied. Moreover, the exterior derivative of the 
fundamental form (sum over repeated indices) 

d$ = ^dE a = I a dE a , (1) 

where I a represent the intensive variables dual to E a , is equivalent to the first law of ther- 
modynamics. 

With only the above information about the equilibrium space £ one can just explore its 
topological properties and the analytic behavior of the fundamental equation $ = Q(E a ). 
In some sense this is what has been done in analytic geometry in the last century. In fact, 
the fundamental equation determines a surface on £ whose geometric properties are related 
to the thermodynamic behavior of the system. For instance, if the thermodynamic potential 
is chosen as the Gibbs potential it is known that the critical points of the surface can be 
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used as a criterion to define phase transitions whose order is directly related to the order 
of the differential derivative of the Gibbs potential at which the critical points exist. The 
study of classical thermodynamics systems is based upon this definition of phase transitions 
and the investigation of sections of the equilibrium space. 

To extract more information from the equilibrium space it is necessary to equip it with 
an additional geometric structure. This can be done, for instance, by introducing a metric 
structure g = g a bdE a dE b on S. Since for a given thermodynamic system there exists only 
one fundamental equation §(E a ) that also determines the properties of S, one would expect 
that the metric g is determined in a unique manner by &(E a ). A particularly interesting 
choice is the Hessian metric 

a* = -**- (2) 
dE a E b ' 

The specific choice with the entropy S as the potential $ can be considered as due originally 
to Rao [7] who in 1945 proposed to use a metric structure in statistical physics, and con- 
structed a particular metric whose components in local coordinates coincide with Fisher's 
information matrix. Rao's original work has been followed up and extended by a number 
of authors (see, e.g., jg] for a review). Later on, in 1975, Weinhold ^ proposed to use 
the internal energy U as the thermodynamic potential. The corresponding metric can be 
shown to be flat for the fundamental equation of the ideal gas, and curved in other more 



general cases, 
by Ruppeiner 



he choice of the entropy as thermodynamic potential was proposed in 1979 
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11] , obtaining a metric which is conformally equivalent to Weinhold's 
metric, with the inverse of the temperature as the conformal factor. The physical meaning 
of Ruppeiner's metric lays in the fluctuation theory of equilibrium thermodynamics. It turns 
out that the second moments of fluctuation are related to the components of Ruppeiner's 
metric. More recently, Liu, Lii, Luo and Shao isj proposed to use the potential $ = U, 
where U is any of the thermodynamic potentials that can be obtained from U by means of a 
Legendre transformation. In the special cases U = U and U = S, one obtains the Weinhold 
and Ruppeiner metrics, respectively. 

The pair (£, g) is then a Riemannian manifold with specific geometric properties that are 
expected to be related to the thermodynamic properties of the system. As mentioned in the 
introduction, it is expected that the curvature of the equilibrium manifold be related to the 
notion of thermodynamic interaction. If this turns out to be true, it is then interesting to 
consider the thermodynamic interpretation of curvature singularities. Since in differential 



geometry the presence of curvature singularities is considered as indicating the break down 
of the geometric formalism, one can expect that such singularities are also related to the 
break down of the thermodynamic description. On the other hand, the limit of applicability 
of classical equilibrium thermodynamics coincide with the occurrence of phase transitions. 
We conclude that curvature singularities should correspond to phase transitions. 

An additional important aspect is the fact that ordinary classical thermodynamics is 



invariant with respect to 
thermodynamic potential 



^egendre transformations, i. e., it is independent of the choice of 
6) . Clearly, the geometric properties of the Hessian metrics defined 
above are not necessarily invariant with respect to Legendre transformations. In fact, since 
a Legendre transformation includes a change of potential $—>•<& that is accompanied by an 
interchange of extensive and intensive variables, it could change the functional dependence 
of the thermodynamic potential $, implying modifications of the explicit components of the 
metric (jSJ). This change of functional dependence could in principle lead to a change of the 
geometric properties of the equilibrium manifold. 

One might then wonder whether the use of non Legendre invariant metric structures in 
ordinary thermodynamics would always lead to results that do not depend on the thermo- 
dynamic potential. Indeed, several examples are known in the literature in which a change 



of thermodynamic potential leads to a modification of the thermodynamic geometry [13Hl8|. 
Some puzzling results and inconsistencies arise also in connection with the use of different 
metrics in the equilibrium manifold 19l-l22j. in the sense that for the same thermodynamic 
system the resulting geometry can be either flat or curved, depending on the choice of 
thermodynamic potential chosen for generating the metric. 

It then follows that the Legendre invariance of g is necessary in order to guarantee in gen- 
eral that the geometric properties of the equilibrium manifold £ do not depend on the choice 
of thermodynamic potential. GTD is a formalism that incorporates Legendre invariance into 
the geometric structure of £ in a consistent manner. Since the Legendre transformations in- 
volve the thermodynamic potential as well as the extensive and the intensive thermodynamic 
variables, to introduce Legendre invariance it is necessary to consider an auxiliary space in 
which all the above variables are independent. This can be done by using the thermody- 
namic phase space T with coordinates Z A = {Q,E a ,I a }, where I a represent the intensive 
variables and A = 0, l,---,2n. Notice that a specific thermodynamic system cannot be 
considered in T, but in the particular subspace where the fundamental equation $ = &(E a ) 
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holds. 

A Legendre transformation can now be defined as a coordinate transformation in T such 
that y 

{Z A } — ► {Z A } = {$, E a , I a } , (3) 

$ = $ - c^ fc /' , = -r, E j = E\ P = E\ V = P , (4) 

where % U j is any disjoint decomposition of the set of indices {1, n}, and k, I = 1, i. In 
particular, for % = {1, ...,n}, i.e., 

$ = $ - 5 ab i? a / b , = -I a , I a = E a , (5) 

we obtain the total Legendre transformation whereas for i = 0, we obtain the identity 
transformation, i.e., 

$ = $, £ a = E a , I a = E a . (6) 

The main point is now that all the geometric structures to be defined on T must be 
invariant under Legendre transformations in such a way that its n— dimensional subspaces 
£ inherit this property. First, we introduce the fundamental 1— form 

= rf$ - 5 ab I a dE b , 5 ab = diag(l, 1, ■ ■ • , 1) , (7) 

which is Legendre invariant in the sense that it transforms as O — > = d$ — 5 a bI a dE h . 
Notice that the existence of this 1— form is guaranteed as a consequence of Darboux theorem 



24J because it satisfies the condition 

6 A (dQ) An j£ , (8) 

which is interpreted as stating that (the tangent space of) T is maximally non-integrable. 
Notice that the (2n + l)-form 6 A (dQ) An = d$ A dl 1 ■ ■ ■ A dl n A dE l ■■■ A dE n determines 
the volume element d 2n+1 Z in T so that the condition (|Sj) indicates that the integral is a 
well-defined operation. The pair (T, 0) is known as a contact space. 

Let us now introduce a Riemannian metric structure G = GABdZ A dZ B in T . From all 
possible metrics G we choose only those that are invariant under Legendre transformations. 
Then the triad (T, 0,G) determines a Riemannian contact manifold that is Legendre in- 
variant and is known as the phase manifold. It is now clear that any geometric structures 
derived from and G under certain conditions can be Legendre invariant. 



Let us now consider the equilibrium manifold £ as a submanifold of T determined by the 
smooth embedding map 

tp : £ T , (9) 

or in coordinates 

99 : {E a } 1 — ► {Z A {E a )} = {&(E a ),E a , I a {E a )} , (10) 

satisfying the condition 

p*(Q) = , i.e., 6 = 5 ab I a dE b on £ , (11) 

which is equivalent to the first law as given in Eq.flTJ). As a result of the smoothness of tp 
(and the pullback p*) and of the invariance of 0, the first law (TTTT) is Legendre invariant in 
this formalism. 

We can now consider the canonically induced metric 

87 A f)7 B 

9 = p*(G)^g ab = —-^G AB , (12) 

as the metric of the equilibrium manifold. In this way, the Riemannian manifold (£, g) can be 
considered as invariant in the sense that it is obtained in a canonical manner from the phase 
manifold (T,Q,G) by using only geometric objects that are Legendre invariant. Notice 
that the definition of the smooth map p as given in Eq. ([TO]) implies that the fundamental 
equation $ = Q(E a ) must be known in order to define the equilibrium manifold. This is in 
accordance with the intuitive description of the equilibrium space given above. 

From the above description it follows that the only freedom in the construction of the 
phase manifold is in the choice of the metric G. Although Legendre invariance implies 
a series of algebraic conditions for the metric components Gab |4j], and it can be shown 
that these conditions are not trivially satisfied, the metric G cannot be fixed uniquely. It 
is important to mention that a straightforward computation shows that the flat metric 
G = 5ABdZ A dZ B is not invariant with respect to the Legendre transformations given in 
Eq. (J4j) . It then follows that the phase manifold is necessarily curved. Moreover, one can 
show that the Hessian metrics for the equilibrium manifold can be generated from a specific 
metric G H of the phase manifold according to 

g H = p*(G H ) = p*(Q 2 + 5 ab dE a dI b ) = g^- b dE a dE b . (13) 
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It is easy to prove that G H is not invariant under Legendre transformations, indicating that 
results obtained by using Hessian metrics could depend on the choice of thermodynamic 
potential. 

If we limit ourselves to the case of total Legendre transformations, we find that there 
exists a class of metrics, 

G = (d$ - I a dE a f + A (UE a I b ) (x cd dE c dI d ) , (14) 

parametrized by the diagonal constant tensors ^ a b and Xab, that are invariant for several 
choices of these free tensors. The constant parameter A is introduced to guarantee that the 
units of the first and second terms coincide. No particular significance is attributed to this 
constant. 

Since the tensors £ a & and Xab must be constant and diagonal, it seems reasonable to ex- 
press them in terms of the usual Euclidean and pseudo-Euclidean metrics S a b = diag(l, 1) 
and rjab = diag(— 1, 1, 1), respectively. Then, for instance, the choice £ a b = 5 a b , Xab = &ab 
corresponds to a Legendre invariant metric which has been used to describe the geomet- 



ric properties of systems with first order phase transitions 



I. 



25]. Moreover, the choice 



iab = S a b, Xab = Vab turned out to describe correctly second order phase transitions espe- 



cially in black hole thermodynamics |25l428l]. At the moment we have no definite explanation 
for the fact that the signature of Xab is able to differentiate between phase transitions of first 
and second orders. Nevertheless, it seems that this difference can be used to propose an al- 
ternative invariant definition of phase transitions. This task is currently under consideration, 
and will be presented elsewhere. 

Obviously, for a given thermodynamic system it is very important to choose the appro- 
priate metric in order to describe correctly the thermodynamic properties in terms of the 
geometric properties obtained in the context of GTD. In this work, we will consider the van 
der Waals gas that is characterized by the presence of first order phase transitions. Then, 
the appropriate metric is 

G=(d<$>- I a dE a f + k(I a E a )(dI h dE b ) , (15) 

from which we can compute the corresponding induced metric 

^•< G > = A (IH ( ' (16) 



It is then clear that once the fundamental equation $(i? a ) is given, the expression for the 
metric of the equilibrium manifold can be found explicitly. 



III. GEOMETRIC PROPERTIES OF THE EQUILIBRIUM MANIFOLD 

Let us consider the fundamental equation for the van der Waals gas in the entropy 
representation 

S = hn(u+^-) +ln{V-b) , (17) 



2 V V, 

where U is the internal energy, V is the volume, and a and b are constants. The parameter b 
is associated with the volume of the molecules of the gas and plays a qualitative role in the 
description whereas a is responsible for the thermodynamic interaction. The corresponding 
intensive variables can be obtained from the first law of thermodynamics dS = ^dU + j;dV 
as 

_ (dS)- 1 _ 2 / a\ (dS\ (dSy 1 _ 2UV 2 - aV + 3ab 

T ~[du) -3{ u+ v)> p -[dv)[du) ~ sv 2 (v - b) ■ (18) 

According to Eq.( fT6|) . in this case the metric 

J dS TT dS^\(d 2 S 1TT2 n d 2 S „ Tnr d 2 S nr2 \ 

" = A [W u + 8v v ) {w du + 2 wsv dUdv + ssn AV ) • (19) 

defines the equilibrium manifold, and for the van der Waals gas reads 



_ A 5 UV 2 - 3 UVb -aV + 3ab 
9 ~ 2 (UV + af(V -b) 

where 



\v 2 dU 2 + 3adUdV-^^dV 2 



(20) 



W(U, V) = 2 V 4 U 2 - 2 V 3 Ua - a 2 V 2 + 12 aV 2 bU + 6 Vba 2 - 6 ab 2 UV - 3 b 2 a 2 . (21) 

In the limiting case a = and 6 = (ideal gas limit) the above expression reduces to the 
metric 

5A /3 dU 2 dV 2 \ 

9 = -Y{2lP + v 2 -) ' (22) 

whose curvature tensor vanishes identically. This proves that the equilibrium manifold of 
the ideal gas is flat, indicating the lack of thermodynamic interaction. 

In the general case a ^ and b ^ 0, the curvature tensor is different from zero. We 
interpret this result as an indication of the presence of thermodynamic interaction. The 
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corresponding curvature scalar can be written as 

«= 5 , ( 23 ) 

(V 3 U -2V 2 a + 6Vba-3 b 2 af (5 UV 2 - 3 UVb - aV + 3 afr) 3 
so that the curvature singularities are determined by the zeros of the two polynomials enter- 
ing the denominator. The function N V( jw(U,V) is a polynomial that is different from zero 
at those points where the denominator vanishes. Using the expression for the pressure given 
in Eq. (fT8l) . one can show that 

V 3 U - 2 V 2 a + 6 Vba - 3 b 2 a =-{V-b) (PV 3 - aV + 2ab) (24) 

and 

5 UV 2 - 3 UVb -aV + 3ab = 3V(V - b)(U + PV) . (25) 

It follows that there exist curvature singularities at those points where the condition PV 3 — 
aV + 2ab = is satisfied. In classical thermodynamics it is known that this condition 
determines the points where first order phase transitions occur in the van der Waals gas 
6j. As for the second polynomial, one can see that no zeros exist for positive values of the 
pressure, a condition that is usually assumed as valid in van der Waals systems. 

We conclude that in GTD the equilibrium manifold of the van der Waals gas correctly 
describes the corresponding thermodynamic behavior. 

IV. GEODESIC INCOMPLETENESS 

In classical thermodynamics, a quasi-static thermodynamic process is a process that 
happens infinitely slowly and therefore it can be ensured that the system passes through a 
sequence of states that are infinitesimally close to equilibrium and, consequently, the system 
remains in quasi-static equilibrium. 

Since each point of the manifold £ represents an equilibrium state, a quasi-static process 
can be interpreted as a sequence of points, i. e., as a curve in S. In particular, the geodesic 
curves of £ along which the laws of thermodynamics are satisfied can represent quasi-static 
processes that take place inside an isolated system without influence from outside. 



If we introduce the concept of thermodynamic length as L = J ds = J \J g a bdE a dE b in £, 

the vanishing of the variation 5L = leads to the geodesic equation 

d 2 E a ^ dE b dE c . . 

r V-7--7- = , (26) 



dr 2 oc dr dr 
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FIG. 1: Geodesies in the equilibrium manifold of the van der Waals gas with initial value V{t = 
0) = 0.1, U(t = 0) = 0, V{t = 0) = 1, and different initial values U(t = 0) in the range [0, 140]. 

where T a bc are the Christoffel symbols of the thermodynamic metric g, and r is an arbitrary 
affine parameter along the geodesic. 

In the case of the equilibrium manifold of the van der Waals gas, the geodesic equations 

are 

d 2 U rr (dU\ 2 rj dUdV ( dV\ 2 . . 



+ r^l^] 2 + 2r^^ + r^ffV = o, (28) 



dT' ' UV \dT J dT dT \dT J 

where the explicit form of the Christoffel symbols can be calculated from Eq. (I20|) . The 
resulting differential equations are highly non-trivial and cannot be treated analytically. 
Instead, we perform a numerical analysis for a large number of different initial conditions. 
The behavior of the geodesies is illustrated in FigJTJ The range of the initial conditions was 
chosen such that the geodesies always reach a maximum value for V(r). 

To clarify this point we perform an analysis in a smaller range V(r = 0) G [0, 14] as 
depicted in Figf2J The main observation is that the geodesies are incomplete, i.e, there 
exists a maximum value of the affine parameter r max for which the numerical integration 
delivers an end value of U max and V max . After this point the geodesies are not defined. To 
understand the origin of this behavior, we analyze numerically the end points U max and V max 
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FIG. 2: Geodesies in the equilibrium manifold of the van der Waals gas for different initial values 
U(t = 0) and V(r = 0) = 0.1. The incompleteness of the geodesies occur at the points where first 
order phase transitions occur. 

and found that they satisfy the relationship 

VL X U max - 2 Vl ax a + 6 V max ba - 3 b 2 a = , (29) 

which, according to Eq.( l24l) . determines the points where the curvature of the equilibrium 
manifold diverges and first order phase transitions take place. It then follows that the 
geodesic incompleteness is due to the appearance of phase transitions. This result cor- 
roborates the fact that curvature singularities can also be defined by using the concept of 
geodesic incompleteness (see, for instance, [291). We conclude that in GTD we can also use 
the geodesies to find the points where phase transitions can occur. 

V. CONCLUSIONS 

In this work, we presented the formalism of GTD starting from the geometric properties 
of the equilibrium manifold and by demanding that they do not depend on the choice of 
thermodynamic potential. In this manner, it becomes clear why it is necessary to introduce 
the auxiliary phase manifold where Legendre transformations are well-defined. Since the 
equilibrium manifold is defined by means of a smooth map as a subspace of the phase 
manifold that is manifestly Legendre invariant, it inherits the invariance property. 

13 



We use a particular metric of the phase manifold, which is is invariant under total Leg- 
endre transformations, to compute the induced metric for the equilibrium manifold. It is 
shown that all the geometric structure of the equilibrium manifold is determined by the 
fundamental equation only. 

The case of a thermodynamic system described by the van der Waals fundamental equa- 
tion is analyzed in detail. First, we found the explicit form of the metric of the equilibrium 
manifold, and showed that it corresponds in general to a curved space, indicating that the 
curvature can be used as a measure of the thermodynamic interaction. In the limiting case 
of an ideal gas, the curvature vanishes as expected for a system with no thermodynamic 
interaction. The curvature singularities of the van der Waals equilibrium manifold were 
shown to be located only at those places where first order phase transitions occur. This, 
of course, is an indication of the break down of the equilibrium thermodynamics approach 
and, consequently, of GTD. 

We investigated numerically the geodesic equations of the van der Waals equilibrium 
manifold. It was shown that for a certain range of initial values, the geodesies cannot be 
extended after a particular value of the afline parameter. The points where this geodesic 
incompleteness was detected turned out to coincide with the points where the curvature 
diverges. This is in accordance with the result known in differential geometry about the 
equivalence between curvature singularities and geodesic incompleteness. We interpret this 
result in GTD as an additional criterion for detecting phase transitions in the equilibrium 
manifold. 
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